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Abstract 

Using the gauge/gravity correspondence, we study the dynamics of a heavy quark in two 
strongly-coupled systems at finite temperature: Super- Yang- Mills in the presence of a mag- 
netic field and non- commutative Super- Yang-Mills. In the former, our results agree qualita- 
tively with the expected behavior from weakly-coupled theories. In the latter, we propose a 
Langevin equation that accounts for the effects of non-commutativity and we find new in- 
teresting features. The equation resembles the structure of Brownian motion in the presence 
of a magnetic field and implies that the fluctuations along non- commutative directions are 
correlated. Moreover, our results show that the viscosity is smaller than the commutative 
case and that the diffusion properties of the quark are unaffected by non-commutativity. Fi- 
nally, we compute the random force autocorrelator and verify that the fluctuation-dissipation 
theorem holds in the presence of non-commutativity. 
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1 Introduction 

Over the years, the study of thermal properties in quantum field theories has been of great 
interest. In the last decade, the discovery of the gauge/gravity correspondence [H El E] has 
provided tools for the study of a large class of strongly-coupled non-abelian gauge theories. 
Aside from pure theoretical motivation, one the the main reasons for pursuing such studies 
is the possible connection with the phenomenology of ultra-relativistic heavy ion collisions 
and/or strongly correlated condensed matter systems. An excellent account of these efforts 
can be found in the reviews [U El El [7] and the references therein. 

The simplest and best-understood example of the correspondence relates four-dimensional 
M = 4 Super- Yang-Mills (SYM) to type-IIB string theory (or supergravity) in asymptotically 
AdSs x S 5 . In particular, the inclusion of a black hole in the bulk is known to be dual 
to a strongly-coupled gauge theory at finite temperature. The prescription to compute 
correlation functions using the gauge/gravity duality was developed a decade ago in [21 E] 
and was subsequently generalized to the real-time finite-temperature formalism in [HE]. This 
tool opened the possibility of computing the decay rates and timescales for the approach to 
thermal equilibrium of certain disturbances [TUl [Til [12] . which in the gravity side are mapped 
to the computation of quasinormal frequencies. Perhaps, the most intriguing result in this 
front has been the well-known result of 1/47T for the shear viscosity to entropy density ratio 
at strong coupling [TSJ, [H] which in turn led to the holographic calculation of many other 
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transport coefficients [151 US] • More recently, it has been possible to gain further insight in 
the previously uncharted out-of-equilibrium regime by considering geometries that evolve in 

time nnnBiiTsiiag. 

In addition, other thermal properties have been inferred by considering various types of 
partonic probes, and analyzing the manner in which the plasma damps their motion. These 
studies include quarks [211 [221 [23], mesons [211 [251 [26], baryons \2T\ [21], gluons [23 [30] 
k-quarks [2Z] and various types of defects [21] • Although all these probes provide different 
information regarding the nature of the plasma, in this paper we focus on the dynamics of 
heavy quarks and their interactions with the thermal bath. 

In the context of this duality, a heavy quark on the boundary theory corresponds to the 
endpoint of an open string that, at finite temperature, stretches between the boundary and 
the black hole horizon. The seminal works [211 [22] focused on the energy loss of a quark 
that is either moving with constant velocity as a result of being pulled by an external force, 
or is unforced but moving nonrelativistically and about to come to rest. Further analyses 
[321 [331 EI] made it clear that this mechanism of energy loss is closely related with the 
appearance of a worldsheet horizon (not to be confused with the spacetime horizon). 

The study of quark fluctuations due to its interaction with the thermal bath involves 
going beyond the classical description of the string. As customary, small perturbations about 
the average embedding are described by free scalar fields propagating on the corresponding 
induced worldsheet geometry. These fields can be excited due to Hawking radiation emitted 
by the worldsheet horizon, which in turn populates the various modes of oscillation of the 
string. It was then found that, once these modes are quantized, the induced motion of 
the string endpoint is correctly described in terms of Brownian motion and its associated 
Langevin equation [351 EE]- This result was obtained by two different approaches: the authors 
of [35] reached this conclusion by assuming (following [37, 38J) that the state of the quantized 
fields is the usual Hartle-Hawking vacuum, which describes the black hole in equilibrium 
with its own thermal radiation. The authors of [36] followed a different but equivalent route, 
employing the relation between the Kruskal extension of the Schwarzschild-AdS geometry 
and the Schwinger-Keldysh formalism [3U1 [9], together with the known connection between 
the latter and the generalized Langevin equation. These calculations were later elaborated 

on in mum mm as Eg. 

In this paper, we generalize the original computation of |35j to the case of non-commu- 
tative Super- Yang-Mills (NCSYM). Non-commutative theories are known to lead to many 
qualitatively new phenomena, both classically and quantum mechanically; in particular, 
the existence of non-local interactions reminiscent of the UV/IR mixing found in string 
theory [16^1 . A simple model of non-commutativity was described by Bigatti and Susskind 
[19] where they considered a pair of opposite charges moving in a strong magnetic field. In 
the limit of large magnetic field, the charges are frozen into the lowest Landau levels and 
the interactions of such particles include the Moyal bracket characteristic of field theories 
on noncommutative space. Similarly, in the context of string theory, it was shown that the 
endpoints of open strings constrained to a D-brane in the presence of a constant Neveu- 

1 For a review of non-commutative quantum field theories, see for example [47l [48] . 
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Schwarz B- satisfy also the non-commutative algebra. A holographic realization of 
NCSYM was given a short time afterwards in [5U1 EE]- 

One of the main motivations to study Brownian motion in the presence of non-commutati- 
vity is the idea that non-local interactions might lead to significant deviations in the behavior 
of the thermal properties of the theory [531 EI]- I n [55] for instance, it was found that the 
rate of decay of a fluctuation propagating in this thermal bath is remarkably larger than 
in the case of ordinary SYM, which leads to faster thermalization. Such a property is 
possibly linked to the connection between the ultraviolet and infrared regimes of the theory, 
which implies in particular that the transverse size of dipoles grows with their longitudinal 
momentum [49J. In order to investigate further properties of this non-commutative system, 
we study in this paper the holographic realization the Brownian motion of a heavy quark. 
More specifically, our aim is to formulate a Langevin equation that accounts for the effects 
of non-commutativity and to study diffusion processes within the plasma. 

The structure of this paper is as follows: in section [21 we describe the gravity dual to 
NCSYM presented in [50| |5T] and, following [35], we review the generalities of Brownian 
motion in the context of the gauge/gravity correspondence. In section [3] we study the 
case of Brownian motion in ordinary SYM with a magnetic field, which is achieved by 
the introduction of a gauge field in the open string sectoio. Apart from being an interesting 
problem in its own rightQ, this exercise helps us to gain some intuition and to set the grounds 
of our computations. In section H] we turn to the study of Brownian motion in NCSYM. We 
begin by postulating a Langevin equation for the non-commutative plasma similarly to the 
one describing the quark in the presence of a magnetic field. We show that this equation 
correctly captures behavior of Brownian particle. This equation is expressed naturally in 
terms of matrices and given its structure, it automatically implies that fluctuations along 
different directions are correlated. We then compute holographically the drag coefficient 
from the response of the Brownian particle to an external force and, finally, we study the 
diffusion process of the quark within the plasma, which turns out to be unaffected by non- 
commutativity. In section El we make some comments about our results and close with 
conclusions. 

2 Background and String Action 

2.1 A gravity dual of NCSYM 

Quantum field theories on noncommutative spaces have been studied intensively in the last 
few years. The essential postulate of non-commutativity is the commutation relation 

[a^z"] =i0"", (1) 

2 By gauge invariance, this is equivalent to a constant magnetic field on the brane. 

3 Here, we should emphasize that we do not expect to recover the non-commutative results in the strong 
B limit. On one hand, there is a critical magnetic field at which Schwinger pair production is energetically 
favored [5 6) . On the other hand, for large magnetic fields the backreaction on the geometry is unsuppressed. 

4 Indeed, many holographic systems present interesting features when a background magnetic field is 
turned on. For a recent review on this topic see [57] . 
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where 6^ is a real and antisymmetric rank-2 tensor. The algebra of functions on a non- 
commutative space can be viewed as an algebra of ordinary functions with the product 
deformed to the non-commutative Moyal product, defined by 



(01*02) 0) 



2 U 



y=z=x 



(2) 



Non-commutative field theories emerge in string theory, as the worldvolume theory of 
D-branes with a constant Neveu-Schwarz 5-field provided that one takes a special limit to 
decouple the open and closed string sectors [581 E3 EHj. Basically, one scales the string 
tension to infinity and the closed string metric to zero while keeping the 5-field fixed. For 
the discussion in the present paper, we will focus specifically on the four-dimensional SU(N) 
NCSYM at finite temperature, taking the non-commutativity parameter to be non-zero only 
in the (x 2 , x 3 )-plane, that is [x 2 ,x 3 ] ~ i9. In the spirit of the gauge/gravity correspondence, 
the dynamics of this theory at large N and at strong 't Hooft coupling should be described 
by a bulk gravitational system. Indeed, a gravity dual for this theory was proposed in [501 EE] 
which, in the string frame, reads 
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u 2 (—f(u)dxl + dx\ + h(u) (dx\ + dxl)) + 



du 2 



u 2 f(u) 



(3) 



u 3 h(u) 



where R 4 = 4ngN = A and g is the value of the string coupling in the IR. In the expressions 
above, 



/(«) = 1-^ 



is the usual Schwarzschild factor and 



h{u) 



1 



1 + a A u A 



(4) 



(5) 



The parameter a which appears in the above expressions is related to the non-commutativity 
parameter 9 of the boundary theory through a = X 1 ^^. Here, A is the 't Hooft coupling of 
the boundary large- iV NCSYM theory. The radial direction u is mapped holographically into 
a energy scale in the gauge theory, in such a way that u — > oo and u — > are respectively the 
UV and IR limits^]. The directions x^ = (t, x) are parallel to the boundary and are directly 
identified with the gauge theory directions. The five-sphere coordinates are associated with 
the global 5*^(4) internal symmetry group, but they will play no role in our discussion. 



5 Throughout this paper, we use the terms "UV" and "IR" with respect to the boundary energy. Then, 
in bulk terms, UV means near the boundary whereas IR means near the horizon. 
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The Hawking temperature of the above solution, which is interpreted as the temperature 
of the non-commutative boundary theory, is given by T = u^j tt. Notice that this temperature 
is the same as the temperature one obtains for the Schwarzschild-AdSs solution, which is 
dual to a thermal state of four- dimensional SYM. Indeed, it is easy to show that all the 
thermodynamic quantities obtained from ([3]) are the same as the ones obtained from the 
Schwarzschild-AdSs solution [51]. In the limit of vanishing temperature (uh — > 0), we are 
left with a geometry that is dual to the vacuum of the NCSYM. The closed string sector 
describing small or large fluctuations on top of this background captures the nonperturbative 
gluonic (+ adjoint scalar and fermionic) physics. The black hole geometry d3J) is a special 
case of such a fluctuation, and is dual to a thermal ensemble of the theory. 

For u <C a -1 , the background ([3]) goes over to the AdSs x S 5 solution. This observation 
just reflects the fact that the non- commutative boundary theory goes over to the ordinary 
commutative SYM theory at length scales much greater than X^-^y/O. On the other hand, 
for u ^> a -1 the background ([3]) shows significant deviation form the AdSs x S 5 solution and 
the bulk spacetime is no longer asymptotically AdS. The boundary theory interpretation of 
this regime just means that the effect of non-commutativity becomes pronounced for length 
scales which are at the order of, or smaller than, A^v^Ij. To summarize, the background 
02]) represents a flow in the boundary theory from a UV fixed point which is NCSYM at 
large N and large A to an IR fixed point given by the ordinary commutative SYM (again, 
at large N and large A). 



2.2 Heavy quarks and holographic Brownian motion 

From the gauge theory perspective, the introduction of an open string sector associated with a 
stack of Nf D7-branes in the geometry ([3]) is equivalent to the addition of Nf hypermultiplets 
in the fundamental representation of the gauge group SU(N), and these are the degrees of 
freedom that we will refer to as quarks. For Nf <C N, we can neglect the backreaction of the 
D7-branes on the geometry; in the gauge theory perspective, this corresponds to working in 
a "quenched" approximation that ignores quark loops. For simplicity, we will take Nf = 1. 
The probe brane covers the four gauge theory directions and is spread along the radial 
direction from u — > oo to u = u m where it ends smoothl)0. 

An isolated quark is dual to an open string that extends radially from the flavor brane 
at u — u m to the horizon at u — The string couples to both the metric and the .B-field 
so its dynamics follows from the action S = Sng + 5*3 

S = [ drdaC = ^— [ drda (j- det g ab + B mn d T X m d a X n ) (6) 
J s 2na' J E V / 



6 Since the reliability of the solution (J3)) requires i? 4 = A to be large, the effect of non-commutativity in 
the boundary theory is visible even at large length scales [51]. This should be compared with the weak- 
coupling regime of the theory where the threshold length scale (beyond which the theory becomes effectively 
commutative) is roughly at the order of \fd. 

7 At this position the S 3 C S 5 that it is wrapped on shrinks down to zero size. 

8 For a review of quark dynamics in the context of the gauge/gravity correspondence see [52] . 

9 The string also couples to the dilaton but this coupling is suppressed by a factor of the string length. 
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where g a b = G mn d a X m di ) X n is the induced metric on the worldsheet E. We choose to 
work in the static gauge, where r = t, a = u. One can easily verify that the embedding 
X m = (t, 0, 0, 0, u) is a trivial solution and this correspond to a quark that is is equilibrium 
in the thermal bath. 

Notice that the string is being described in first-quantized language and, as long as it is 
sufficiently heavy, we are allowed to treat it semiclassically. In gauge theory language, we 
are coupling a first-quantized quark to the SYM fields, and then carrying out the full path 
integral over the strongly-coupled fields, treating the path integral over the quark trajectory 
in a saddle-point approximation. The mass of the quark m is related to the position of the 
flavor brane u m and can be obtained by a straightforward computation: 



m 



2-rca' 



du y/g tt g uu = ^-{u r 

271 



. R 2 

u h) ~ TT U ^ 
2ix 



for u m > u h . 



(7) 



Now, we want to study fluctuations over the above embedding. In particular, we are 
interested in motion on the Moyal plane, so our ansatz for the perturbations will be X m = 
(t,0,X 2 (t,u),X 3 (t,u),u). Using this, the induced metric has the following components: 
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Quu 

9tu 
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a'R 2 u 2 



-f + h[X 2 2 +Xl 
h [X 12 + X' 2 ) + 



w 4 / 



a'R 2 u 2 h [X 2 X' 2 + X 3 X' 3 ) , 



(9) 
(10) 



where Xi = dtXi and X[ = d u Xi. Up to quadratic order in the perturbations, the action 
can be written as 



S 



I? 
Air 



J dtdu u 4 fh [X' 2 + X'i) - j [±l + X 2 ) + 2a 2 u 4 h (x 2 X' 3 - X 3 X' 2 



'ir. 



Note that we dropped the constant term that does not depend on Xi. We can consider 
also the situation in which one has a forced motion due to a electromagnetic field in the 
background. This can be easily realized by turning on a world-volume U(l) gauge field on 
the flavor brane. Since the endpoint of the string is charged, this amounts to add the minimal 
coupling to the action S = Sng + <Sb + <Sem, where 



5 EM = / (A t + AiX t )dt, 

JdE V J 



'121 



This will exert the desired force on our heavy quark. However, this coupling is just a 
boundary term, so it will not play any role for the string dynamics in the bulk, other than 
modify the boundary condition. We shall ignore this part of the action for now but we will 
come back to it later on. 

Because t is an isometry of the background (j3]), we can set 



X{t,u) ~ e-^gM 



(13) 
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and use the frequency u to label the basis of solutions to the equations of motion. Since the 
action (II ip is quadratic in the perturbations, we expect linear differential equations. The 
solutions are particularly easy to obtain near the horizon limit u ~ Uh, where the action 
reduces to 

S « u 2 h h(u h )^ J dtdu* \{X'i + X' 2 ) - [xl + X 2 ) + 2a 2 u 2 h (x 2 X' 3 - X 3 xf) 

w u 2 h(u h )^ J dtdu* [(X> 2 2 + X' 2 ) - (Xl + X 3 2 )] . (14) 

Here, the primes denote derivatives with respect to the tortoise coordinate u*, which is 
defined by 

du * = 277 v ( 15 ) 

Note also that the last term drops out because it is a total derivative. Thus, the equations 
of motion are then 

K - df)X t = 0, (16) 

which show that in this region Xi behave like massless Klein-Gordon scalars in flat space. 
The two independent solutions are 

X, (out) ( M ) = e- iut gl ont \u) ~ e --(*-^) (17) 

and 

X f n \u) = e-^gf n \u) ~ e -^ t+u *\ (18) 
corresponding to outgoing and ingoing waves respectively. Near the horizon one finds that 

M l,og(^") (19) 



up to an additive numerical constant, so 

, (out/in) f s _ / U 



gV nt ' m >(u) ~ 1^- - lj . (20) 

Away from the horizon, these solutions will have a complicate dependence, but it still holds 
that c/(° ut ) = g^* (see appendix lAl). 

Standard quantization of quantum fields in curved spacetime [65J leads to a mode expan- 
sion of the form ^ 

Xi(t,u)= / — [a u u u (t, u) + aluu(t,u)*\, (21) 



2tt 1 

where the functions u u (x) correspond to a basis with positive-frequency modes. These modes 
can be expressed as a linear combination of outgoing and ingoing waves with arbitrary 
coefficients, i.e., 

Uufau) = A [g (ont) (u) + Bg^ a \u)] e~ lujt . (22) 
7 



The constant B is fixed through the boundary condition at u = u m but one generally 
obtains that it is a pure phase B = e ld (see sections I3T31 and 133]) . The outgoing and ingoing 
modes have, then, the same amplitude and this implies that the black hole, which emit 
Hawking radiation, can be in thermal equilibrium [68]. The constant A on the other hand, 
is obtained by requiring the normalization of the modes through the standard Klein-Gordon 
inner product. 

For any functions fi(t, u) and gj(t, u) satisfying the equations of motion, the Klein-Gordon 
inner product is defined by 

where a is a Cauchy surface in the (t, u) part of the metric, g is the induced metric on a 
and n M is the future-pointing unit normal to a. It can be shown that this inner product is 
independent of the choice of cr |65j, but for simplicity we take it as a constant-t surface. 

We want to normalize using (I23j) . However the main contribution to the integral 
comes from the IR region [42J, which in terms of the tortoise coordinate is just 

R 2 f 

(fi,9j)a = -iSij—ulh(u h ) du*(fig* - fog*). (24) 

Of course, there is a contribution to the inner product from regions away from the hori- 
zon, but because the near- horizon region is semi- infinite in the tortoise coordinate u*, the 
normalization of solutions is completely determined by this region. 
After some algebra, we find that 



A -]/coRV h h(u h y (25) 

so that (m w ,m w ) = 1, ensuring that the canonical commutation relations are satisfied: 

[a w , Ou>] = [al, a*,] = 0, [a u , a*,] = 2n5(u - u'). (26) 

In the semiclassical approximation, the string modes are thermally excited by Hawking 
radiation emitted by the worldsheet horizon. In particular, they satisfy the Bose-Einstein 
distribution: 

,t „ A = 2^-0/) 



Using this and the mode expansion given in (|2ip . we can derive a general formula for the 
displacement squared of the Brownian particle. 

First of all, let us identify the position of the heavy quark as the string endpoint at the 
boundary u = u m , i.e., 

x t (t) = Xi(t,u m ) = [ ^ [a^u^t^Um) + aluufaum)*] . (28) 
Jo Zlx 
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Then, it follows that 



f°° dudu' r 

(xi(t)xi(0)) = / (a a; a^)M tJ (t,M m )M a; /(0,M m )* + (a^a aJ /)M aJ (t,M m )*M tJ /(0,« r . 

Jo ( 27r ) L 

(29) 

This has an IR divergence that comes from the zero point energy, which exists even at zero 
temperature. To avoid this we simply regularize it by implementing the normal ordering 
: a u al : 



la^a^:, and after doing so we ge 
(:Xi(t)xi(0):) = [ ^ e p}_ 1 [v>w(t, u m )u u {Q, u m )* + u u (t, u m )*u u (0, u m )\ 



MW£*W \ g ^(u m ) +Bg i-Hu m )\\ (30) 
From this correlator, we can compute the displacement squared of the quark as: 

where we have replaced the explicit dependence on A. For future reference we also compute 
the general form of the momentum correlator, 

(:p ( (f)Pi(0) : ) = -m 2 8?(:n(i)n(0) : ), 

d^AA^ut) |g(out) K) + Bs(1 „ )(Mm)r (32) 



3 Brownian Motion in SYM with a Magnetic Field 

Non- commutative SYM theory has a dual interpretation in terms ordinary SYM with a 
large and constant magnetic field. We start by studying the Brownian motion in this second 
system which, aside from being an interesting problem in its own right, helps us to gain some 
intuition and to set the physical grounds of our computations. Although we find that there 
are some similarities between these two configurations, our final results show that there are 
some features that are qualitatively different. 



3.1 Langevin dynamics in the presence of a magnetic field 

The problem of the Brownian motion of a charged particle in an external magnetic field was 
first investigated almost fifty years ago in the seminal papers [HH E2] ■ This is an old topic that 
has originated a lot of interest and it is of great importance in the description of diffusion and 
transport of plasmas and heavy ions. Nowadays, together with the free Brownian motion, 
it is widely used as a classic textbook example of how transport properties and correlation 
functions should be computed in generic situations governed by the Langevin equation. 

10 Another way to regularize it is to use the canonical correlator introduced as in 35 . However, this does 
not change the late-time or low- frequency behavior of the correlator. 
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The discussion of this section will be around the field theory description of Brownian 
motion in the presence of a magnetic field. Later on, we will show how to realize this 
phenomenon at strong-coupling, in terms of a probe string living in a black hole background. 

Let us consider the Langevin equation of a charged particle of mass m and unit charge 
q, in presence of a magnetic field B: 

P(t) = -7oP(t)+v(t) xB + R(t), (33) 

where p(i) = mv(t) is the momentum of the Brownian particle and v(t) = x(£) its velocity. 
The terms in the right-hand side of fl33|) correspond to the friction, lorentz force and random 
force, respectively, and 7 Q is a constant called the friction coefficient. One can think of the 
particle as moving under the influence of the magnetic field, losing energy to the medium 
and at the same time, getting random kicks as modeled by the random force. 

As a first approximation, we can assume that the random force is white noise, with the 
following averages: 

(Ri(t)) = 0, (RityRjtf)) = K 5 i3 5(t - £'), (34) 

and where k q is a constant which, due to the fluctuation-dissipation theorem, is related to 
the friction coefficient through 

7o = (35) 

This is due to the fact that the frictional and random forces have the same origin at the 
microscopic level, i.e., collisions with the particles of the thermal bath. 

If the magnetic field B = Bx is pointing along the x-direction, we can write (1331) in the 
matrix form 

p i (t) = -A ijPj (t) + R i (t), (36) 

where 

j lo \ 
Kj = 7o -u . (37) 
\ w„ To / 

Here u = B/m denotes the Larmor frequency. Since the magnetic field is oriented along 
the x-direction, it affects the motion in the y and z directions only. Fluctuations along 
the x-direction decouple and are unaffected by the presence of the magnetic field. We thus 
restrict our attention to the fluctuations in the transverse plane. 

In order to decouple the remaining equations we have to diagonalize the the above matrix. 
The normal modes are p± = P2 ± ip3 with corresponding eigenvalues A± = 7 =F iw . Thus, 
defining x± = £2 ± ^3 and R± — Ri ± iRz we get 

P±(t) = -\±p ± (t) + R ± (t), (38) 

whose formal solution is given by 

p± (t) = e" A± V±(0) + f e- x±( t-^R±(f)dl!, (39) 

Jo 
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x±(t) - x ± (0) = — - (l-e- A± > ± (0) 
TJIA4 



f R±(t')dt' - f e- x[t - tn >R±{t')dt' 
Jo Jo 



(40) 



Using the above, we can immediately obtain x 2 (t) , xs(t) by taking the real and imaginary 
parts of x±(t). p2{t) and pz(f) can also be obtained in a similar way. In thermal equilibrium, 
the two-point correlation functions of P2{t) and ps(t) are given by 

(p 2 (t)p 2 (0)) = ^e-^cos(u; t), (41) 

(p 3 (*)i> 3 (0)) = (p 2 (t)p 2 (t')) , (42) 
(P2(t)p 3 (0)} = ^e-^smiuot). (43) 

The relevant point here is that, unlike in the case of free Brownian motion, now the com- 
ponents p 2 (t) and ps(i) are correlated. But, the autocorrelator (pi(t) 2 ) of each individual 
component has the same value ~ n /2 , -f = mT as in the case of zero magnetic field. This 
can be easily understood by recognizing that (pi(t) 2 ) is twice the kinetic energy multiplied by 
the mass of the particle and that this quantity does not change by application of a magnetic 
field. Thus, the time scale associated to the energy loss due to drag, t re iax ~ V7°> * s then 
independent of the magnetic field. 

Another important scale related to the Brownian motion of the quark is the time asso- 
ciated to diffusion. This can be derived by computing the two-point correlation function of 
x 2 (t) and xs(t), from which one can infer the following late-time behavior for the displace- 
ment squared: 

s 2 (t) = {[ Xi (t) - ^(O)] 2 ) ~ 2Dt, for t»l/ 7o . (44) 

where D is called the diffusion constant. In the presence of a magnetic field, the fluctuation- 
dissipation theorem leads to 

D = O 2 2 1 2 = — t 2> ( 45 ) 

which decreases with increasing magnetic field. Thus, the relevant time scale for correlations 
~ D is smaller in this case and this implies that the diffusion process is more efficient. 

The Langevin equation, as presented in (133]) . captures the essential properties of the 
stochastic processes in Brownian motion, but it fails to give a physically consistent picture 
for sufficiently short times t, in which the particle suffers only a few or no impacts. It is a 
general feature of any dynamical system that the dynamical coherence becomes predominant 
in short time scales, or at high frequencies. Thus, we are led to a natural extension of the 
Langevin equation in the form [631 El] 



P(t) = - f df7(*-*')p(0+v(t) xB + R(f) + E(*), (46) 
J —00 

where 

(Ri(t)) = 0, (RiWRjffl) = mAt - f )■ (47) 
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The main difference between the generalized Langevin equation (|46|) and the usual one (1331) 
is that the friction depends now on the past history of the particle through 7(t), called the 
memory kernel, and that the random forces at different times are not independent. Note 
that we have also included a fluctuating external force E(t) that can be applied to the system 
(e.g., an external electric field). 

For a magnetic field pointing in the x-direction B = Bx and focusing on the transverse 
fluctuations, we get 

P±{t) = - [ dt'^it - t>±(0 ± ^oP±it) + R±{t) + E±(t), (48) 

J — oo 

where p± = p 2 ± ip 3 , R± — Ri ± iR 3 , E± = E 2 ± iE 3 and u Q = B/m. 
In frequency domain, the above equation is simplyM 

7[wJ — too =F 

and taking the statistical average of the same we get 

(P±M) = lJ-±{v)E±(u), with /x±(w) = -p-j 7 — • (50) 

The quantity fi±{u>) is known as the admittance. Thus, we can then determine the admittance 
fj,±(cu), and thereby j[cj], by measuring the response (p±(cu)) to an external fluctuating force. 
In particular, if the external force is taken to be 

E ± {t) = e~^K ±1 (51) 

then (p±(t)) is just 

ip±(t))=fji ± {u)e- iut K ± = fM ± (u)E ± (t). (52) 

For late times (or low frequencies) the generalized Langevin equation reduces to its non- 
local progenitor, and the timescales associated to the decay of the two-point function of 
the momentum as well as the displacement squared are the same as discussed before. In 
particular, note that 

/i ± (0) = l — (53) 



11 Causality imposes that j(t) — for t < so j[uj] in this expression denotes the Fourier-Laplace transform, 

poo 

7 H=/ dtj^e^, 
Jo 

while p±(u>), R±(uj), and E±(uj) are Fourier transforms, e.g., 

/•oo 

j»±(«)= / dt P± (t)e^. 
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so 



7 ° = ^ = Re Uo)J' (54) 

For a quantity O, the power spectrum Io{u) is defined as 

/■oo 

Io(uj)= / dt(O(t )O(t + t))e^\ (55) 



and it is related to the two-point function through the Wiener-Khintchine theorem 

(O(u)O(u')) = 2tc5{u + u/)I (u>). (56) 

For stationary systems ( 151)j) is independent of t , so one can set t = in such situations. 
When the external force is set to zero, from (I49D it follows that 



p ± (co) = R±{UJ) , = K")R±(to), (57) 
7[cuJ — zoo =F «w„ 

and, using ( 156]) one gets that 

I p± (u) = \^u)\ 2 I R± (u). (58) 
Therefore, the random force correlator appearing in (1471) can be evaluated as 

« ± ( w ) = /^( w ) = -^j. (59) 

This will be important in the next section, to check the validity of the fluctuation-dissipation 
theorem. 

3.2 Bulk dynamics and the drag coefficient 

We now turn to the holographic realization of Brownian motion in the presence of a magnetic 
field. Let us begin by considering the action (II ip. which in the commutative limit a — > 
reduces to 

?2 r r . . _ i / . . . _n i 

(60) 



Sk^ Jdtdu u'f{X'i + X'i)-- f {xl+Xl 



This action describes the dynamics of a string in Schwarzschild-AdSs and it is dual to a 
quark in ordinary SYM at finite temperature. We then turn on a gauge field in the flavor 
brane of the form 

-> B 

A=—(yz-zy) } (61) 

thus getting the desired magnetic field B = Bx. As explained before, this only appears as a 
boundary term (Tl2|) . so it will not affect the bulk dynamics of the string. The equations of 
motion coming from the above action are: 

= fd u (u 4 fXl) -Xi. (62) 
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We now proceed by expanding X; in modes as in (fT3l) . i.e., 

Xi(t,u) = e-^giiu). 
Then the equations of motion ( |62l) can be written as 



r-i 

where we defined dimensionless quantities 

u 

Uh' 



y 



(y 4 



;9i(y) 



(63) 



(64) 



(65) 



and where primes denote now derivatives with respect to y. The wave equation for the 
modes (1641) is independent of the magnetic field and it is exactly the same as the equation 
considered in [35] for d = 4. 

We need to find the solutions of the equation ( J64l) . In general, it is not possible to 
do this analytically for an arbitrary frequencies v and hence we employ a low frequency 
approximation C 1 by means of the so-called matching technique. Here, we only write 
down the final result, relegating the details of the computation to appendix [K\ The two 
solution that correspond to outgoing and ingoing waves at the horizon behave asymptotically 
as 
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(out/in) 



(y) 



1 =F -u(tt - log(4)) 1 + _ T _ + 0(l/y 



2y< 



w 



(66) 



We now consider the forced motion of our Brownian particle due to a fixed external 
magnetic field and a fluctuating electric field. As mentioned in section 12.21 this amounts 
to the addition of the boundary term fTl~2l to the action ([6]), which imposes a boundary 
condition of the form 



« \dE 



OX' 



(67) 



BE 



Here, = —(F it + FijX^) is the usual Lorentz force and 



( 

-E x 
-E 2 
\-E 3 





-B 2 



E 2 
-B 3 




E3 
B 2 
-B x 




(68) 



/ 



Then, for a magnetic field pointing in the x-direction, B x = B, and transverse electric fields 
E 2 = E 2 (t) and E 3 = E 3 (t) we find a set of boundary conditions that inevitably mix the 
fluctuations along the transverse directions: 



d2 



BX, 



Eo 



^u'fX', + BX 2 



E 3 , 



(69) 



where u m denotes the position of the flavor brane. Our goal is now to compute the thermal 
expectation value (or one-point function) of the momentum, and then extract the admittance. 
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The general solution for X{ is the sum of ingoing and outgoing modes at the horizon Xi = 
^(out) x(out) + ^(in) X (in) 5 where X (out/in) = e -iwt g (out/in) _ In the sem i c i ass i ca i approximation, 

outgoing modes are always thermally excited because of Hawking radiation, while the ingoing 
modes can be arbitrary. However, because the radiation is random, the phase of v4 J - out ' ) takes 
random values and, on average (^4j out ^) = 0. Then, we can write (Xj) = (A\ )e~ ujt g^(u), 
where correspond to the normalized solution given by (16"6"j) . For the remaining part 

of this section we will denote (Af n ^) = Ai and = g for simplicity. 

In the Brownian motion literature, it is customary to work in a circular basis when an 
external magnetic field is included. Thus, we define 



X± = X 2 ± iX 3 = e- lU]t g±(u) and E± = E 2 ± iE 3 = e~ luJt K±. 
In fact, the equations (jSHD decouple in this basis. After some algebra we get 

R 2 



2tt 



u 4 fg'A ± ±uBgA ± 



K 4 



where A± = A 2 ± iA 3 . Inverting this relation we obtain 

K ± 



A, 



%u*fg'±uBg 



(70) 



(71) 



(72) 



from which we can read the average position of the heavy quark, (x±(t)) = (X±(t,u m )) 



-iujt 



g(u m )A±, 



(s±(t)) 



-iuit 



§u*fg'±uBg 



and 



<P±(*)> = E±(t) 



—tumg 



^u*fg'±uBg 



This expression allows us to read the admittance, 

— iujmg 



H±\u) 



^u*fg'±LoBg u=Um 

In the zero frequency limit, and for a heavy quark u m ««(,, we obtain 

2m 



//±(0) 



iry/XT 2 ± 2iB' 



from which we can infer 



nVXT 2 , , B 

7o = — and u = ±— . 

2m m 



(73) 
(74) 

(75) 

(76) 
(77) 



As expected, the friction coefficient is not modified by the presence of the magnetic field, 
which is consistent with the fact that the magnetic field does not do work. Also, we find 
that 



^relax 

7o 



2m 



(7E 



15 



This temporal scale dominates the late-time decay of the one-point function of p(£), for 
a quark that transverses the plasma, in agreement with the previous works |2TJ [22]. The 
late-time^ behavior is dominated by the low frequency limit of the generalized Langevin 
equation, case in which f H6|) reduces to its nonlocal progenitor ( }33l) . From (1391) we can thus 
infer that, 

(p±(t)> ~ e - 7oi e ± ^'. (79) 

This is exactly what is expected for the Brownian motion of a charged particle in the presence 
of a magnetic field. To obtain the thermal averages of P2(t) and Ps(t) we can simply take 
the real and imaginary parts of (1791 . 



3.3 Diffusion and the fluctuation-dissipation theorem 

The purpose of this section is to compute holographically the displacement squared of the 
heavy quark and to extract from it the diffusion constant D. The upshot of the computation 
is summarized in equation (13TI) . which is valid for an arbitrary background. However, the 
functions fu (u) as well as the details for the computation of the constant B vary according 
to each situation. 

Before proceeding with the direct calculation of this quantity, it is useful to understand 
the boundary conditions we want to impose on the fields. Although we are interested in 
the world-sheet theory of the probe string, in the static gauge the induced metric on the 
string inherits the geometric characteristics of the spacetime background. This means that 
the usual rules for correlators in the gauge/gravity correspondence apply in our case. 

First, we need to impose an UV cutoff in order to have a quark with finite mass. The 
natural place to impose the cutoff is given by the location of the flavor brane u m , which 
can be related to the mass m of the quark through (J7J). The mass of the quark is chosen to 
be dominant scale of the system, so usually one would push the cutoff up to the boundary 
u m — > oo and choose normalizable boundary conditions for the modes. However, in our case 
that would correspond to a infinitely massive quark and there would be no Brownian motion. 
A Neumann boundary condition at u — u m also does not work in our case because we would 
go back to the case of free Brownian motion. Instead, we use a mixed boundary condition 
in which the external magnetic field is on but the fluctuating electric fields are turned off. 
According to (1B"9~1) this is, 



I? 
2tt 



u A fX' 2 - BX 3 



0, ^fX' 3 + BX 2 
In 



0. 



(80) 



or in terms of the modes X±, 



^u 4 fX' ± ±tBX ± 



0. 



51) 



12 There should be a smooth crossover with the early-time regime or high frequency limit. See for example 



67 
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Recall that X± can be expressed as the sum of outgoing modes and ingoing modes found 
previously, with arbitrary coefficients. Following the convention introduced in (122 p . let us 
write 

X±(t,u) = At [g {mt) (u) +B±g<M(uj\ e**. (82) 
From ( 18 ip it follows then that 



B+ = 



= e ie± . (83) 



E?_ u 4f g (in), ± Wj B^(in) 

The fact that B± is a pure phase is self-evident, given that 

^(out) = ^(in)* 

. To leading order 

in frequency one finds that 

B± = ;w r -|iT^ + om = + od/^ + oh, (84) 



7rR 2 T 2 (y? n -l)±2iBy4 l v y \irR 2 T 2 ±2iB 
from which one gets 

| g (o„ )+B±9Mr = _^g_ + OM . (85) 

The late-time behavior of the displacement squared can be then inferred from the low- 
frequency limit of ( 13~TT) . i.e., 



16\/AT 3 r 00 , sin 2 (f) 4ttVAT 3 

Thus, as expected, we find that the diffusion constant defined as in is given by 

n = 2?rv/XT3 = T -Tq (z 7 ) 
4B 2 + vr 2 AT 4 m 7 2 + w 2 ' ^ 

Finally, in order to give an explicit check of the fluctuation dissipation theorem (135)) we 
compute the random force autocorrelation appearing in ([59]) to extract the coefficient k . 
From ([32]) we can evaluate the two-point correlator of the momentum p± as follows: 

du 2m 2 u 2 \A\ 2 cos{ut) , (out) (in) ,2 



/ \ /„s v ^ 2m to \A\ cos(cot) , /■ mit \ / N „ r;ni 
:ft=(*)A=(0):> = I ^ e ^ 1 > W mt \u m ) + B ± g^ 



f 00 —— ^L\n^(u ) + B + a^(u )Fe** 



Therefore, 



^w=(w) = ^ T e /J|^ i b ( ° ut) («m) +B ± gW(u m )\\ (89) 

and combining this with ([75]) one finds that 

J R± ( W ) = ttv^T 3 + (90) 

This gives us precisely a coefficient k = ir\/XT 3 which agrees with ( ]35l) . providing and 
explicit check of the fluctuation-dissipation theorem in the presence of a magnetic field. 
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4 Brownian Motion in NCSYM 



We now turn our attention to the study of Brownian motion in the non-commutative setup. 
The main difference here is that the closed string sector is modified by the inclusion of an 
antisymmetric £>-field. In this setup, and after the appropriate decoupling limit, the effective 
field theory is described by a gauge theory living in a noncommutative space. It is interesting 
to study the similarities and differences with our previous computation when the magnetic 
field was introduced in the open string sector. 

4.1 Langevin dynamics in the non- commutative plasma 

To start, let us postulate a generalized Langevin equation of a particle in a non-commutative 
thermal bath: 

Piit) = ~ [ dt'Yij{t - t') Pj {t') + Riit) + Ei{t), (91) 

J — oo 

where 

(Ri{t)) = 0, (Ri^Rjit')) = Kij {t - t'). (92) 

In this case, the 5-field does not appear explicitly in the Langevin equation, though its 
effect should be somehow be present through the coefficients V and k. We propose that in 
this case T is a matrix that encodes the effects of the non-commutativity. In particular, for 
non-commutativity in the (x 2 , x 3 )-plane, we propose to write 

( \ 

T tJ (t) = j(t) -n(t) . (93) 

In the low frequency limit, the above equation becomes local in time, allowing us to write 

Pi (t) = -Tijpj (t)+Ri(t) + Ei(t), (94) 

where ^{t - t') = ^y 5(t - t'), Q(t - t') = tt 5(t - if) and n{t - if) = K 8{t - If). Note 
that this has exactly the same structure of (|36j) . with 7 G being the usual friction coefficient 
and fl playing the role of the Larmor frequency. Furthermore, if the fluctuation-dissipation 
theorem applies, we expect that the relation ( 13 5 p is also true for the present configuration. 
Solutions fl39l) and (j40p hold in the low frequency limit. This means that the two-point 
correlators fj4Tj) - fj43|) . as well as the diffusive behavior of the displacement squared (jSJ), act 
in the exact same way, but now with new coefficients j Q , Q D and n D which might depend on 
the noncommutative parameter 9. 

The fluctuations along the x 1 -direction are unaffected by the presence of the non-commu- 
tativity. We thus restrict our attention to the fluctuations on the Moyal plane. These 
fluctuations can be decoupled by working in the circular basis p± = P2 ±^P3, R± = Ri ±iRj, 
and E± = E 2 ± iE s . The eigenvalues of ( 193|) are found to be X± = 7 =F iQ, so we can rewrite 
(ED as 

P±{t) = - [ dt' A±(t - t')p ± {t') + R±{t) + E±{t). (95) 
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In frequency domain this equation can be written as 



R ± (u) + E±(u) 
A ± [oj\ — lu 

and taking the statistical average we obtain 

(P±M) = fi±{u)E±(uj), with n±(w) = —r^ — — . (97) 

Then, by measuring the response (p±(u))) due to an external force we can determine the 
admittance fi±(u)) and thereby /i±[u;]. In particular, if the external force is taken to be 

E ± (t) = e~ iu>t K ±) (98) 

then 

(p±(t)} = fi±(uj)e- i " t K± = ti±(u)E ± (t). (99) 

With this at hand, one can take the real and imaginary parts of fi± to obtain 7 and Q 
respectively. Note also that in the zero frequency limit 

MO) = y-fnT = -— ^- (100) 

The analysis of the power spectrum and two-point functions is the same as the one done in 
section 13.11 In particular, the equation ( 159]) relating the random force correlations with the 
momentum correlations should apply in this case, which will be useful to check the validity 
of the fluctuation-dissipation theorem for the current setup. 

4.2 Bulk dynamics and the drag coefficient 

For the action given by (11 ip we can derive the following equations of motion: 

= Ld u (u A fhX' 2 - a 2 u 4 /iX 3 ) - d t (X 2 - a 2 u 4 fX^ , (101) 
= Ld u (u A fhX! 3 + a 2 u A hX 2 ^j - d t (A 3 + a 2 u 4 fX!^ . (102) 
The term with mixed derivatives cancels out in both cases and one ends up with 

= i-d u (u 4 fhX' 2 ) - Aa 2 u 3 fhX 3 - X 2 , (103) 
f 

= J -d u {u 4 fhX' 3 ) + 4a 2 u 3 fhX 2 - X 3 . (104) 

These are two coupled partial differential equations. We now proceed by expanding X i in 
modes by setting 

X 2 (t,u) = e-^g 2 (u), and X 3 (t, u) = e^~^g z {u). (105) 
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We have introduced a phase difference for reasons that will become clear below. The equa- 
tions of motion for the modes are 



= ^d u (u 4 fhg' 2 ) +Atua 2 u 3 fhg 3 e^ + u 2 g 2 , (106) 
= ^d u (u i fhg' 3 )-4iua 2 u 3 fhg 2 e-^ + u 2 g 3 . (107) 

If we choose e tLp = ±z, or equivalently ip = ±7r/2, the two equations of motion turn out to be 
the same. This motivates to consider the linear combinations X± = X 2 ± iX 3 = e~ luJt g±(u), 
where 

9+ + 9- , 9+ ~ 9- nna s 

92 = and g 3 = — — — . (108) 

Not surprisingly, this is completely equivalent with the circular basis introduced in section 
|3] for the case of Brownian motion in a magnetic field. In this basis the equations of motion 
decouple, and can be rewritten as: 

= ^ + ( ^-i)(i + 6V) ^ (?/) + \W^W 2 (^-i)(i + ^v) J " ±(y) (109) 

where we have defined dimensionless quantities 

u UJ 

y = — , v = — , b = au h , (110) 

and the primes denote derivatives with respect to y. The normal modes X± correspond to 
fluctuations with circular polarization, rotating clockwise or counterclockwise, respectively. 

Explicit solutions to the above equations can be found in appendix IA.21 The final ex- 
pressions for outgoing and ingoing modes are 



(in) 



and 



# >to> ~ U - + (l + + W - log(4))) (l - £ ) + 0(l/v* 



b 2 ±i\ 2y 2 J V h2±l 8 VV % 

(112) 

respectively. 

— * 

We now exert an external fluctuating force E(t) on the string endpoint by turning on an 
electric field Fu = E\ on the flavor brane. Variation of the whole action implies the standard 
dynamics for all interior points of the string, but now with boundary condition 



nrL^H = £„ (113) 



where E; is the external force. 
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From ( ITTj) . it follows that 



— (u 4 fhX' 2 - a 2 u 4 hX 3 ) = E 2 , — (u 4 fhX' 3 + a 2 u A hX 2 

2w \ J 2ir \ 



(114) 



where u m denotes the position of the D7-brane. Our goal is to find the admittance of the 
system for which we need the one-point function of the momentum. Again, it is convenient 
to work in the circular basis, i.e. X± = X2 ± iX^ and E± = E2 ± 1E3. The general solution 
for X± is the sum of outgoing and ingoing modes X± = ^fa ut )jQ° ut ) _|_ A^X^. However, 
as discussed before, the phase of A± nt ^ takes random values and on average (4<° ut >) = 0. 
Then, we can write (X±) = (A± )e~ Mt g±(u) and E± = e~ wt K± but, for simplicity, we will 



denote (A±^) = A± and g K ±' = g± in the remaining part of this section. In this basis the 
boundary conditions decouple: 

R 2 



2tt 



(u 4 fhA±g' ± (u) ± ua 2 u 4 hA ±g± (u)) 



K 



±- 



;ii5) 



or, in terms of the dimensionless quantities y, v and b as defined previously, 

2ttK± 1 + 6V 



A. 



y=y m 



R 2 u\ (y*-l)g> ± (y)±vb 2 y*g ± (y) 
The average position of the heavy quark, (x±(t)) = (X±(t,y m )) is then 

2tt (i + bY)g±(y) 



(116) 



and 



(x ± (t)) 



(P±(t)) 



R 2 u\ {y 4 -l)g' ± {y)±vb 2 y 4 g ± {y) 
—2irium (1 + b 4 y A )g±(y) 



K±e 



Uijt 



y=Vm 



from which we can read the admittance, 

—27riiom 

H±{u) - 



(y 4 -l)g' ± (y)±ub 2 y 4 g ± (y) 

(i + bY)g±(y) 



y=Vn 



R 2 u\ (y 4 - l)g' ± {y) ± vb 2 y±g ± {y) 



(117) 



;ii8) 



(119) 



y=yv 



In the zero frequency limit and for large massif! we get 



J"±(0) 



2m(l + b 4 yf n ) \^ib 2 



2m 



b 2 RHT 2 b 2 y^ ± i R 2 ttT 2 



[lTib 2 ) + 0{l/yi 



(120) 



As we can see, the real part coincides with the expected value for the commutative case, but 
now there is an additional part that is imaginary (and independent of the temperature, given 



13 In the large mass expansion, we set the non-commutative parameter to a fixed value and then take the 
limit y m — > 00. 
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that b = auh = airT). For a quark that transverses the plasma one can read the following 
evolution at late times: 

(P±(*)) ~ e-^e ±in °*, (121) 

where 

nVXT 2 7T 3 AgT 4 
7o ~ 2m(l + 7r 4 A# 2 T 4 ) and ""Ml+^W)' 

The friction coefficient is modified by the presence of the non-commutativity, and in this 

case 

_ 1 _ 2m(l + 7 r 4 Ag 2 T 4 ) 
relax "7o" v^T 2 ' ( j 



This agrees with a previous computation of the drag force in the Maldacena-Russo back- 
ground [69l [70] (in the non-relativistic regime) and implies that the non-commutative plasma 
is less viscous in comparison to the commutative one. 



4.3 Diffusion and the fluctuation-dissipation theorem 



Now we turn to the computation of the displacement squared in the Maldacena-Russo back- 
ground. First of all, we need to understand the boundary conditions we want to impose on 
the fields. In this case, the effect of non-commutativity is already present in the background 
itself, so the free Brownian motion is realized by imposing a Neumann boundary condition 
at u = u m . According to (I114p . this means 



or simply 



u^fhX' 2 - a 2 u A hX : 



0. 



2tt 



fX'± ± ia 2 X± 



u 4 fhX' 3 + a 2 u A hX 2 



0. 



0. 



(124) 



(125) 



where again, we defined X± = X2 ± iX^. The general solution is then written as a linear 
combination of outgoing and ingoing modes, 



X ± (t,u)=A ± [g^ t \u) + B ± g^\u) 
and from the boundary condition (I125P we get that 



-iujt 



(126) 



±ua 2 g { l nt) 



(in) 



J9± 



Here, B± is also a pure phase given that, in this case, it still holds that g± ut ^ = g±^ 
leading order in frequency, we obtain 

{l T ib 2 ){l±ib 2 y i m 



B, 



(l±ib 2 )(lTib 2 yi 
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(127) 
\ To 

(128) 



With this at hand, we can also compute 



t} + B ± g^\ 2 = (j^±^ + 0(l/y* )J + C{v\ (129) 

and finally, by taking the low-frequency limit of f l3Tj) we compute the late-time behavior of 
the displacement squared: 

„ , N 16 f°° 7 sin 2 4 . . 

Note that the factor that depends on b exactly cancels with the h(uh) term that appears 
in the normalization constant A. We find that, surprisingly, the diffusion constant is not 
affected by non-commutativity, but the relation with respect to 7 D and Q Q is still the same 

2 T ^ 

D = ^^ = . (131) 

This suggests that the fluctuation-dissipation theorem holds even in the presence of non- 
commutativity. In order to explicitly check this, we compute the random force correlator in 
to extract the coefficient k q . From (1321 it follows that 



(:p±{t)p±m = £ t T^T ] e£\ - 1 \9 (ont \n m )+B ±9 ^(u m )\ 2 e-^ (132) 
and hence 

= m ^J ] i \g^\nm) + B ±g W(u m )\\ (133) 

At leading order, we find through f l5^|) that 



ttVXT 3 
l + vr 4 Afl 2 T 4 ' 



(134) 



This agrees with ( )35|) . thus providing an explicit check of the fluctuation-dissipation theorem 
for the non-commutative plasma. 



5 Discussion 

In this paper we carried out an analytical study of the dynamics of a heavy quark in two 
strongly-coupled systems at finite temperature: SYM in the presence of a magnetic field 
and NCSYM. The former was realized by studying the fluctuations of a string living in an 
AdS black hole background and turning on a gauge field in the open string sector. The 
latter was achieved by replacing the background for one that incorporates the effects of non- 
commutativity through the introduction of an antisymmetric .B-field in the closed string 
sector. 
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For both systems, we found that the Langevin equation that describes the dynamics of 
such a quark has matrix coefficients and this fact induces correlations along the relevant 
directions. This is in complete agreement with the classical theory of Brownian motion in a 
magnetic field [6H [62]. We then displayed the basic properties of these equations by com- 
puting holographically the admittance and the random force autocorrelator and we showed 
that these two quantities are related through the usual fluctuation-dissipation theorem. The 
existence of such theorem is due to the fact that, at the microscopic level, friction and 
random forces have the same origin, i.e., interactions with the degrees of freedom of the 
thermal bath. Finally, we studied the diffusion of the quark in both systems and we showed 
that, although the non-commutative plasma is less viscous, the late-time behavior of the 
displacement squared is unaffected by the non-commutativity. 

As explained in the introduction, one of the main motivations that led to this work was 
to establish whether the fast thermalization found in [55] holds for more general situations or 
not. An important difference between the approach of this paper and that of |55j is that here 
we studied the non-commutative plasma with a local probe. In the previous work, on the 
other hand, we considered composite non-local operators that are obtained by smearing the 
usual gauge covariant operators over open Wilson lineal The fast thermalization and large 
decay rates of the modes are then possibly related to the non-local character of the probea 15 !. 
It would be interesting to further explore this question by probing the theory with probes of 
different 'size' and study the associated timescales for the approach to thermal equilibrium. 
Two interesting possibilities to consider are Wilson loops and entanglement entropy. 

It is important to emphasize that all of our computations were performed in the low 
frequency limit of the theory, in which case the analytical computations were under control. 
Going beyond the hydrodynamical regime might also offer new insights but it requires a 
numerical approach. For example, in [71] it was shown that a large class of holographic 
quantum liquids exhibit novel collective excitations that appear due to the presence of a 
magnetic field. At high frequency, the dominant peak in the spectral function is associated 
to sound mode similar to the zero sound mode in the collisionless regime of a Landau Fermi 
liquid. The study of Brownian motion within this regime is beyond the scope of this paper, 
but it is also left for future works. 

In conclusion, the results obtained in this paper shed additional light on the thermal 
nature of non-commutative gauge theories and suggest future directions of research. The 
present study constitutes yet another illustration of the usefulness of the gauge/gravity 
correspondence. 
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A Solutions for the string embedding 

In this appendix we well derive explicitly the solutions to the equations of motion considered 
in sections 13.21 and 14.21 At low frequencies, the solutions can be obtained by means of the 
matching technique [351112] (see also [66J ) . To find these solutions, consider three regimes: 
(A) the near horizon solution (y ~ 1) for arbitrary u, (B) the solution for arbitrary y but 
!/< 1, and (C) the asymptotic solution (y — > oo) for arbitrary v. The idea is to find the 
approximate solutions for each of the three regimes, and to match these to leading order in 
v. 



A.l Solution for the AdS-Schwarzschild black hole 

Here we will solve the equation (|64p . by considering the three regimes alluded above. 
(A) In this regime we can focus on the equation as y — > 1: 

o = g"(y) + ^jg'iv) + ^^)- (135) 

We have dropped the subindex i because the equations of motion are the same for both 
i — 2, 3. The general solution in this regime is 



ff 



\ y ) = A(° nt \y - l)-/ 4 + A^ n \y - l)"*/ 4 , (136) 



where the coefficients A^ out ^ and correspond to outgoing and and ingoing modes respec- 
tively. Normalizing each solution according to (120]) and expanding for low frequencies, we 
obtain 



9 



= (y- 1)^ /4 ~ 1 ± j log(y - 1) + 0(v 2 ). (137) 

(B) For the regime of low frequencies we proceed to expand the solution as a series of 
the form 

g B (y) = 9o(y) + m(y) + ^2(2/) + ... (138) 

The first function can be obtained analytically by solving the equation 

= gZ(y) + ^- [ g' (y). (139) 

The general solution goes as follows: 

g (y) = B 1 + B 2 (tan-%) + tanh" 1 ^)) , (140) 

where, in order to have a reliable expansion in frequencies, we have to assume that the 
constants B 1 and B 2 are independent of v. In order to find the appropriate coefficients Bi 
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and B 2 to obtain the ingoing and outgoing modes we expand around the horizon and match 



with the first term in (11371) . After doing so we obtain B\ 

gt t/in) (y) = i. 



(out /in) 



1 and B, 



(out/in) 



0, so 



(141) 



The equation for gi(y) turns out to be the same as for go(y), but now the matching has to 
be done with the second term in (11371) . At the end we get 



(out /in) . 



(142) 



9i 1 \v) = ±^((2 + z)7r + 2zlog(2)) T ^(tan- 1 ( Z /)+tanh- 1 ( Z /)), 

~ ±^i°g(j/ - 1) as y !■ 

Then, up to this order we can write 

^(out/in )(y) = ! ± 1 ((2 + t)n + 2l log(2)) T * (tan- 1 ^/) + tanh- 1 ^)) + G(z/ 2 ). (143) 



Asymptotically, these solutions behave as 



9 



_B(out/in) 



(v) 



1 T -^(tt - log(4)) + (9(z/ 2 



+ 0(l/y±). (144) 



(C) The general solution in region C can be found perturbatively, as a series expansion 
in 1/y. The leading order terms go as follows: 



9°{y) = c x 



1 + v + ° (1/!,4) 



9i 

y 3 



10y 2 



0(l/y 4 



(145) 



Again, comparing this with (11441) . in the low frequency limit, we obtain 



(out/in) _ 1 I ! „„j ^(out/in) ' 



Izp -^(7r-log(4)), and C^ 



=F-z/. 
3 



(146) 



Thus, the normalized asymptotic solutions for modes corresponding to outgoing and ingoing 
waves at the horizon, can be written as 

9° {mt/in \y) ~ (l T - Iog(4))) (l + ^) =F ^ + 0(1/Z/ 4 )- (147) 
This agrees with the solutions reported in [35]. Note in particular that g^ ont '(y) = g^ n '{y)*- 



A. 2 Solution for the Maldacena-Russo background 



Next, we solve (11091) following a similar procedure as the one used for the commutative case. 
(A) In this regime we can focus on the equation as y — > 1: 



i , . , f 2 



;9±(y), 



(148) 
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which is equivalent to (I135p . As expected, the IR is not affected by the non-commutativity. 
The general solution in this regime is 

gi{y) = At at \y - l) w/4 + A^\y - 1)"*/* (149) 

where the coefficients A± and correspond to outgoing and ingoing modes respectively. 
Normalizing these solution and expanding for low frequencies, we obtain 

gi (out/in \y) = (y- i) ±W4 ~ 1 ± j tog(y - l) + 0{y\ (150) 

(B) In this regime we start by expanding the solution as a series in v: 

g±{y) = 9o±{y) + vgi±{y) + v 2 92±{y) + ... (151) 

The first function can be obtained analytically by solving the equation 

4(l + 6 4 )y 3 „ , . , 

o = %±\y) + (y4 _ 1)(1 + 6V) W(y)- (152) 

The general solution goes as follows: 

9o±(y) = B ±1 + B ±2 [b'y - 1(1 + b A ) (tan- 1 ^) + tanh" 1 ^))] , (153) 

where, in order to have a reliable expansion in frequencies, we have to assume that the 
constants B±i and B± 2 are independent of v. We then expand near the horizon and match 
the above solution with the first term in (|150p to obtain the outgoing and ingoing modes. 
We find that 5i7 t/in) = 1 and Si° 2 ut/in) = 0, so 

9 ^ t/in \y) = I- (154) 

Plugging (11511) into the equation of motion and using (11541) we can can derive the following 
equation for the next leading order term of the normalized modes: 

4(1 + b A ) y 3 Ab 2 y 3 
° = 9l± " iy) + (y*-l)(l + b^) 9l± ' iy) ± (y 4 -i)(i + 6V) ' (155) 

The solution to this equation is 

9i±(y) = B ±1 + B ±2 [b A y - |(1 + 6 4 ) (tan- 1 ^) + tanh" 1 ^))] (156) 

(tan- 1 ^) +tanh- 1 (y)) , 



and after expanding around the horizon and matching with the second term in (11501) we 
obtain 



D( out) _ J±z 2 ^ n ± ^ 2 + - ^ 8 - lQ g( 4 ))) + lQ g( 4 ) n M) 

> V±2 ~ h 2 {h 2 Tl y ( i5 ^ 
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and 



R (in) _ , (1 ~ 2i)n Tb 2 {{2 + i)n - 1(8 - log(4))) + log(4) (in) _ i 
±x ~ 8(b 2 ±i) ' ±2 ~ b 2 (b 2 ±i) 



(15* 



Asymptotically, the solutions g± 



B (out /in) , 



(v) = 9iT /in \y) + vg\T /m >(y) + 0[v 2 ) behave as 



(out/in) 



B(out)/ x 



ib 2 v 
b 2 =f z 



y+ 1 



1 



6 2 =f« 8 



- -ii/(7r - log(4)) + £>(z/) 



and 



B(in), x 

g± (y) 



ib 2 v 



:y + l + 



ib 2 v 
b 2 ±i 



+ -iv{ix - log(4)) + C(z/ : 



+ 0(1/1/), (159) 



+ 0{l/y). (160) 



(C) The general solution in region C can be found as a series expansion in 1/y. The 
leading order terms are: 



g±(y) = c± iy 



v 



i-^ + o(iV) 



±2 



(161) 



Comparing this with ( I159P and (I160p we obtain 

ib 2 v 



(out) 



/^r(out) 

ft 2 ?? G±2 



* 2// 1 • / 1 ,A\\ 

l-^-8-(--log(4)), 



and 



(in) 
±1 



CS } = 1 + 7^7 + ^-^(4))- 



(162) 
(163) 



b 2 ±i' ±z b 2 ±i 
Thus, the normalized asymptotic solutions, in the low frequency limit, and for modes corre- 
sponding to outgoing and ingoing waves at the horizon, can be written as 



C(out) 

9± 



(y) 



and 

C(in)/ x 

g± \y) 



ib 2 vy f 


v 2 


b 2 Ti \ 


2y 2 


ib 2 vy ( 


V 2 


b 2 ±i\ 


" V 


9t ut \y) = 


(in) 

9± 



1 - 



%b 2 v 



b 2 Ti 8 



- -iufr - log(4)) 1 - — + 0(l/y 



6y< 



(164) 



ib 2 v 



b 2 ±i 8 



v 



+ -u/(tt - log(4)) 1-FT )+0(W). 



6y s 



(165) 



It is clear that g± (y) = g± (y)* as was found for the case analyzed in appendix lA.il 
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